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Abstract 

X 

$— i ' A comparison of experimental heavy residue cross sections from the 

reactions 86 Kr+ 64 Ni, 112 ' 124 Sn with the model of deep-inelastic trans- 
fer ( DIT ) is carried out. A modified expression for nucleon transfer 
probabilities is used at non-overlapping projectile-target configura- 
tions, introducing a dependence on isospin asymmetry at the nuclear 
periphery. The experimental yields of neutron-rich nuclei close to the 
projectile are reproduced better and the trend deviating from the bulk 
isospin equilibration is explained. For the neutron-rich products fur- 
ther from the projectile, originating from hot quasiprojectiles, the sta- 
tistical multifragmentation model reproduces the mass distributions 
better than the model of sequential binary decay. In the reaction with 
proton-rich target 112 Sn the nucleon exchange appears to depend on 
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isospin asymmetry of nuclear periphery only when the surface separa- 
tion is larger than 0.8 fm due to the stronger Coulomb interaction at 
more compact di-nuclear configuration. 

Introduction 

Nucleus-nucleus collisions in the Fermi energy domain exhibit a large variety 
of contributing reaction mechanisms and reaction products ( see e.g. PQ, [2], 
PI, |3] ) and offer the principal possibility to produce mid-heavy to heavy 
neutron-rich nuclei in very peripheral collisions. In the reactions of massive 
heavy ions such as 86 Kr+ 124 Sn [5 j and 124 Sn+ 124 Sn jH], an enhancement was 
observed over the yields expected in cold fragmentation which is at present 
the method of choice to produce neutron-rich nuclei. Further enhancement of 
yields of n-rich nuclei was observed in the reaction 86 Kr+ 64 Ni in the very 
peripheral collisions, thus pointing to the possible importance of neutron and 
proton density profiles at the projectile and target surfaces. 

In this article, we present an investigation of the reaction mechanism of 
the very peripheral nucleus-nucleus collisions with an emphasis on the pos- 
sible role of different proton and neutron densities at the surface as a factor 
influencing the process of nucleon exchange. When nuclei with different neu- 
tron to proton ratio collide peripherally and nucleon exchange takes place, 
one observes equilibration of the N/Z-ratio. An isoscaling study was carried 
out on the reactions 86 Kr+ 112 Sn and 86 Kr+ 124 Sn at 25 AMeV |3] and a cor- 
relation of the isoscaling parameter with isospin equilibration was observed 
[7j. Peripheral nucleus- nucleus collisions can be described theoretically using 
the model of deep-inelastic transfer [H] in combination with an appropriate 
model of de-excitation. Very good description of experimental observables 
was obtained [21 EJ El HI] using the DIT model of Tassan-Got [TU] and the 
de-excitation codes SMM [TT] and GEMINI |I2j. SMM represents the statis- 
tical model of multifragmentation ( SMM ) and GEMINI invokes the model 
of sequential binary decay ( SBD ). An enhancement of the yields of n-rich 
nuclei over the prediction of such calculations was observed in the reaction 
86 Kr+ 64 Ni at 25 AMeV j2j. The shapes of the velocity spectra suggested a 
process with a short timescale such as the very peripheral collisions where 
the details of neutron and proton density profiles at the projectile and target 
surfaces can play a significant role. An attempt was made in [5] to implement 
the density-dependent correction into the DIT model and an improvement 
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of agreement with experimental data from the reactions 86 Kr+ 64 Ni, 112,124 Sn 
was observed for some regions, which however was not consistent enough to 
be conclusive and a need for further development of a self-consistent density- 
dependent DIT model or an equivalent approach was inferred. 

Model analysis of heavy residue data 
Reaction 86 Kr+ 64 Ni 
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Figure 1: Experimental mass distributions ( symbols ) of elements with Z = 30 - 35 
observed in the reaction 86 Kr+ 64 Ni at 25 AMeV [2] compared to the results of the standard 
DIT calculations combined with the de-excitation codes GEMINI and SMM ( solid and 
dashed line, respectively ). 



In Fig. n we present the experimental mass distributions of elements 
with Z = 30 - 35 observed within the separator acceptance in the reaction 
86 Kr+ 64 Ni at 25 AMeV compared to the results of the standard DIT cal- 
culations [TU] combined with two de-excitation codes SMM jTTj ( dashed line 
) and GEMINI [12J ( full line ), used for de-excitation of the quasiprojectiles 
emerging after the DIT stage. The simulated yields were corrected for angu- 
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lar acceptance of the separator positioned at 0° ( covering polar angles 1.0 
- 2.7 ° [2J ). One observes an excess of experimental yields of neutron- rich 
nuclei with Z = 34 - 32, which is evident both for the SMM and the GEMINI, 
in the former case the effect appears to extend to lower atomic numbers. As 
suggested in [2], such excessive yields of neutron-rich nuclei can be caused 
by the effect of the neutron-rich surface of the target nucleus, which in pe- 
ripheral collisions can lead to stronger flow of neutrons from the target to 
the projectile ( or flow of protons in the opposite direction ), thus reverting 
the flow toward isospin equilibration. 




-10 -5 5 10 -10 -5 5 10 
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Figure 2: a) Calculated dependence of R n — R p on the difference of neutron and proton 
chemical potentials \i n — fx p , obtained using the extended Thomas- Fermi model JS]- Open 
and solid symbols represent the nuclei with A<100 and A>100, respectively, b) As in a), 
but the values of R n — R P are calculated using the liquid-droplet parametrization IT3] . 

In order to understand how the properties of nuclear surface can revert the 
isospin flow one first needs to clarify whether and to what extent such effect of 
nuclear surface is implemented in the standard DIT model of Tassan-Got [TU] . 
The model of Tassan-Got is based on various macroscopic parametrizations. 
For an estimate of the thickness of neutron-skin a parametrization based on 
liquid-droplet model of ref. [T3] is used. The values of neutron-skin thickness 
t are then used to determine the radius of the sharp uniformly charged sphere 
which is used to describe the single-particle Coulomb potential Further- 
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more, the radius of the nuclear part of single-particle potential is additionally 
enhanced by 0.5 fm in order to match the typical radius of the proximity po- 
tential [Tl|, which allows to describe the reaction dynamics properly. Thus 
such adjustement essentially leads to further enhancement of the thickness 
of neutron-rich surface. However, the relation between neutron and proton 
single-particle potentials and corresponding density distributions, on which 
the nuclear surface is usually defined, is not simple and can be, in principle, 
determined only in self-consistent manner. In any case, the DIT model of 
Tassan-Got [TU] proved successful in describing the bulk N/Z-equilibration 
and thus explanation to the enhancement observed in Fig. ^will be appar- 
ently related to the effects not considered in it, specifically these which will 
lead to modification of isospin-asymmetry at nuclear periphery. Also it is of 
interest to identify observables characterizing such effects. 

The observable which is of particular interest is the isospin asymmetry at 
the nuclear periphery, specifically its deviation from the bulk isospin asym- 
metry of the nucleus. It can be shown that the value of this observable at 
the mean half-density radius (R n + R p )/2 is determined dominantly by its 
linear correlation to the the thickness of a neutron-rich surface R n — R p . 
An estimate of the thickness of a neutron- rich surface R n — R p , where R n , 
R p are the half-density radii of neutrons and protons, can be obtained us- 
ing nuclear structure calculations. A calculated dependence of R n — R p 
on the difference of neutron and proton chemical potentials \x n — /i p , ob- 
tained using the extended Thomas- Fermi code of Kolomietz et al. [To] , 
is presented in Fig. |2l along with the liquid-droplet estimate R n — R p = 
2.714((1 - 2Z/A) - 0.0016582Z/A 1 / 3 )/(l + 3.45/A 1 / 3 ) [TU1 [T3|. The calcu- 
lations were performed for the /3-stable nuclei with masses A = 16 - 238 ( 
representing typical target nuclei ) with the neutron- rich isotopes 94 Kr, 132 Sn 
being added in order to extend the N/Z-range. For the extended Thomas- 
Fermi calculation, one observes a linear correlation for both subsets of nuclei, 
for heavier nuclei ( solid symbols ) the values of R n — R p appear to be some- 
what larger than for lighter nuclei ( open symbols ), thus implying a thicker 
neutron-rich surface for heavier /3-stable nuclei. The values obtained using 
the liquid-droplet parametrization ( squares ) follow similar trend, however 
the absolute values are significantly larger than in extended Thomas-Fermi 
calculation. The spread of the values is caused by the use of fi n — \i v from the 
previous plot, which do not follow the trend of isospin asymmetry so strictly 
as it is assumed in the liquid-drop model, where both R n — R p and /i n — fi p 
follow the isospin asymmetry linearly. The global correlation of R n — R p with 
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ggests the possibility to relate the effect of isospin asymmetry at 
the surface to the difference of chemical potentials. However, it should be 
noted that the values of fJ, n —fJ, p calculated using the extended Thomas- Fermi 
model represent the "smooth" ( macroscopic ) nuclear properties and may 
differ from the experimental values of the difference of neutron and proton 
separation energies —(S n — S p ) which is typically used as the estimate of 
— Hp for real nuclei at ground state. On the other hand, when assum- 
ing that the correlation of R n — R p with //„ — \i v is preserved also for real 
nuclei, the experimental values of —(S n — S p ) can be assumed as represent- 
ing the values of R n — R p of such nuclei, with the density profiles modified 
due to the microscopic structure, including the effects not considered in the 
liquid-droplet parametrization used by DIT model of Tassan-Got such as 
shell structure. 
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Figure 3: Experimental mass distributions ( symbols ) of elements with Z = 30 - 35 
observed in the reaction 86 Kr+ 64 Ni at 25 AMeV [2] compared to the results of the modified 
DIT calculations combined with the de-excitation codes GEMINI and SMM ( solid and 
dashed line, respectively ). 



The correlation of the R n — R p with fi n — fi p shown in Fig. El suggests 
the possibility to implement a corresponding correction into the DIT model 
taking into account the effect of nuclear periphery. Based on the discussion 
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of Fig. |21 in the previous paragraph, in particular considering the possibility 
to estimate the surface properties of real nuclei using the value of —(S n — S p ), 
we made a minor modification in the DIT code of Tassan-Got ^U] by scaling 
the transfer probabilities by the exponential factors 



P n (P — > T) > e -o.5K(6s n p-5Spp-6s n T+fiSp T ) p n (p — > 

p ^p y p^ y e 0.5K(SS n p-8Spp-5S n T+SS P T) p ^p > p^ 

p (p y p^ y e -0.5K(8S nT -SSpT-SS n p+SSpp) p (rp ^ p\ 

p (p y p^ y Q0.5K(SSn.T—SSpT—SS n p+SS p p) p (p y p^ 

where k is a free parameter ( due to essentially a Boltzman factor struc- 
ture, possibly representing an inverse statistical temperature ) and 5S n p, 
5S p p, 5S n T, SSpT represent the differences of neutron and proton separation 
energies for the projectile and target calculated using the experimental 16 
and liquid-drop [T7] masses, thus expressing the effect of the microscopic 
structure. The smooth part is subtracted from the experimental values due 
to the fact that the macroscopic values of \i n — \i p follow the bulk N/Z-ratios 
of reaction partners and the bulk N/Z equilibration is described consistently 
by the DIT code of Tassan-Got. Thus, the DIT model is supplemented with 
phenomenological information on shell structure at the nuclear periphery 
which can explain the deviation of nucleon exchange from the path toward 
isospin equilibration, and the model framework assumes the structure com- 
monly used e.g. to describe the ground state properties of nuclei and fission 
|18j . where the liquid-drop model ( or extended Thomas- Fermi model in more 
recent works as for instance in ref. ^J] ) is supplemented by additional term 
describing the effect of shell structure. In the present case, such a term is 
phenomenologically related to shell corrections in neutron and proton sep- 
aration energies SS n p, 5S p p, 5S n T, SS p t, which are based on experimental 
information and thus their use can result in enhanced predictive power. The 
modified DIT calculation was used only for non-overlapping projectile-target 
configurations, consistent with the assumption that it represents an effect of 
nuclear periphery. A cut-off was set at zero half-density surface separation 
( representing touching half-density surfaces ) below which a standard DIT 
calculation, following the path toward isospin equilibration, was used. 

In Fig. El we present the experimental mass distributions of elements 
with Z = 30 - 35 observed within the separator acceptance in the reaction 
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Figure 4: Experimental mass distributions ( symbols ) of elements with Z = 30 - 35 
observed in the reaction 86 Kr+ 124 Sn at 25 AMeV |S] compared to the results of the 
standard DIT calculations with GEMINI ( full line ) and SMM ( dashed line ). 



86 Kr+ 64 Ni at 25 AMeV compared to the results of the modified DIT 
calculations, again combined with the two de-excitation codes GEMINI ^2] 
( full line ) and SMM ^T] ( dashed line ). The simulated yields were fil- 
tered for angular acceptance as in the case of Fig. Several calculations 
were performed with different values of k and the value of k = 0.53, used 
in the modified DIT calculations presented in Fig. El was obtained as an 
optimum value best reproducing the experimental mass distributions. One 
can observe that the modification of the DIT code allows to dramatically 
enhance the agreement with the experimental yields of neutron-rich nuclei 
with Z = 35 - 32 when using both de-excitation codes. The GEMINI cal- 
culation results in the nearly symmetric mass distributions which appear to 
overestimate the widths of mass distributions of lighter elements. The SMM 
calculation appears to reproduce well the yields of neutron-rich nuclei also for 
lighter elements, on the other hand the yields of /3-stable isotopes appear to 
be overestimated. Such differences of GEMINI and SMM calculations are in 
good agreement with the results of the work jl] where the SMM calculations 
lead to better reproduction of yields originating from the hot quasiprojec- 
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tiles, while for the colder quasiprojectiles with excitation energies 1-2 AMeV 
GEMINI performed better due to the implementation of sequential binary 
decay which is missing in the simulation of secondary emission from the hot 
fragments in SMM. 
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Figure 5: Experimental mass distributions ( symbols ) of elements with Z = 30 - 35 
observed in the reaction 86 Kr+ 124 Sn at 25 AMeV compared to the results of the 
modified DIT calculations with GEMINI ( full line ) and SMM ( dashed line ). 



Reactions 86 Kr+ 112 124 Sn 

In order to further examine the modification of the DIT code, the results of 
simulations were compared to data from other reactions. Figs. 01 GO show the 
experimental mass distributions of elements with Z = 30 - 35 observed within 
the separator acceptance in the reaction 86 Kr+ 124 Sn at 25 AMeV jHj again 
compared to the results of the standard ( Fig. |U) and modified ( Fig. EJ) DIT 
calculations combined with the GEMINI ( full line ) and SMM ( dashed line 
) . The simulated yields were filtered for angular acceptance of the separator 
positioned at 4° ( covering polar angles 2.7 - 5.4 ° [5J ) with appropriate 
azimuthal corrections. One can see that in this case, the overall agreement 
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with the experimental data is improved in the modified DIT calculation ( 
using k = 0.53 as in the previous case ) when comparing to the standard one. 
As in the previous case, the modified DIT calculation with SMM appears to 
reproduce the shapes of mass distributions more consistently ( except the 
overestimation of the yields of /3-stable nuclei ), while the GEMINI code 
appears to lead to more symmetric mass distributions with overestimated 
width at lower atomic numbers. 
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Figure 6: Experimental mass distributions ( symbols ) of elements with Z = 30 - 35 
observed in the reaction 86 Kr+ 112 Sn at 25 AMeV [H] compared to the results of the 
standard DIT calculations with GEMINI ( full line ) and SMM ( dashed line ). 



Figs. El present the experimental mass distributions of elements with Z 
= 30-35 observed within the separator acceptance in the reaction 86 Kr+ 112 Sn 
at 25 AMeV again compared to the results of the simulations, the stan- 
dard ( Fig. IH1 ) and modified ( Fig. ) DIT calculations combined with the 
GEMINI ( full line ) and SMM ( dashed line ). The simulated yields were fil- 
tered for angular acceptance of the separator positioned at 4° ( as in the case 
of reaction 86 Kr+ 124 Sn ). In this case the modified DIT calculation ( using the 
value k = 0.53 successful in previous cases ) combined with GEMINI leads to 
improvement for Z = 35 - 34, consistent with previous cases. However, both 
the modified and, to a lesser extent, the standard DIT calculation combined 
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Figure 7: Experimental mass distributions ( symbols ) of elements with Z = 30 - 35 
observed in the reaction 86 Kr+ 112 Sn at 25 AMeV |S] compared to the results of the 
modified DIT calculations with GEMINI ( full line ) and SMM ( dashed line ). Dash- 
dotted line represents the modified DIT with a cutoff shifted to 0.8 fm combined with 
SMM. 



with GEMINI appear to overestimate the yields of neutron-rich nuclei with Z 
= 33- 30 due to shifted centroids and overestimation of both maximum value 
and width. The modified DIT calculation combined with SMM appears to 
reproduce the shapes of mass distribution consistent to previous cases, well 
at the neutron-rich side ( specifically for the Z = 35 - 34 ), while overesti- 
mating the yields of /^-stable nuclei. For Z = 33 - 31 the overall elemental 
yields appear to be overestimated ( even when subtracting the excess yields 
of /3-stable nuclei ) and the standard DIT calculation combined with SMM 
reproduces these yields better. Such a situation in the case of the reaction 
86 Kr+ 112 Sn may signal that an increase of proton transfer probability into 
such a massive proton-rich target according to Eqn. (JTJ) may be reverted at 
smaller separation distances by the effect of increasingly repulsive Coulomb 
interaction. The dash-dotted line in Fig. [7| represents a modified DIT cal- 
culation ( k = 0.53 ) where the cut-off is shifted to minimal separation of 
half-density surfaces equal to 0.8 fm, again combined with SMM. Such a cal- 
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culation reproduces the experimental mass distributions much better, thus 
suggesting that for the proton-rich target 112 Sn, the proton transfer barrier 
assumes its sensitivity to isospin asymmetry of nuclear periphery 0.8 fm out- 
wards when compared to neutron-rich targets, due to the effect of stronger 
Coulomb repulsion at more compact di-nuclear configurations than predicted 
by the approximation used. 

Conclusions 

In summary, the DIT model of Tassan-Got ^H] is supplemented with a phe- 
nomenological correction introducing the effect of shell structure on nuclear 
periphery. A consistent agreement with experimental data is achieved in the 
reactions of a 25 AMeV 86 Kr beam with three different target nuclei, specif- 
ically allowing to describe the deviation of the nucleon exchange from the 
path toward isospin equilibration. The value of parameter k appears to be 
system-independent, consistent with the system-independent correlation of 
neutron skin thickness with difference of neutron and proton chemical po- 
tentials. The success of the modified DIT calculation can be explained as a 
correction reflecting the modification of neutron and proton transfer prob- 
abilities, most probably of the transfer barriers, in peripheral collisions due 
to the effect of shell structure on isospin asymmetry at nuclear periphery. 
Such an effect is of interest for better prediction of the production rates of 
exotic nuclei with a wide range of N/Z-ratios at the new generation of rare 
isotope beam facilities. Discrepancies observed for the most proton-rich tar- 
get nucleus 112 Sn signal a loss of sensitivity of the transfer barriers toward 
isospin asymmetry at the nuclear periphery due to increased Coulomb re- 
pulsion at non-overlapping di-nuclear configurations with surface separation 
below 0.8 fm. The de-excitation code GEMINI appears to systematically 
lead to overly symmetric mass distributions with increasingly overestimated 
widths toward the lighter elements. The SMM appears to reproduce the 
experimental shapes better ( except for excess yields of /3-stable isotopes ). 
In particular the trends of the yields of neutron-rich nuclei are reproduced 
consistently. 

This work was supported through grant of Slovak Scientific Grant Agency 
VEGA-2/5098/25 and by the Department of Energy through grant No. DE- 
FG03-93ER40773. 
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